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1. Introduction 

A precise measurement of gauge boson production cross-sections for pp scattering will be 
crucial at the LHC. W and Z bosons will be produced copiously, and a careful measurement 
of their production cross-sections will be important in testing the Standard Model more 
rigorously than ever before, and uncovering signs of new physics which may appear through 
radiative corrections. In addition, these cross-sections have been proposed as a "standard 
candle" for measuring the luminosity through a comparison of the measured rates to the 
best theoretical calculations of the cross-section. Investigation of this means of measuring 
luminosity began at the Tevatron and will continue at the LHC [jl], Qj. 

In a previous paper, || three of the authors analyzed the systematic errors in Z 
production using state-of-the-art computational tools. A comparable analysis of systematic 
uncertainties in W production appeared in Ref. {|. Since that time, both the experimental 
approach to the measurements and the theoretical results needed to calculate them have 
both been refined, so we will extend a similar analysis as in Ref. || to the case of the W. As 
in the case of the Z, NNLO QCD calculations of these processes, previously available only 
for the total cross-section |3] and rapidity distribution ||, are now available in differential 
form permitting an analysis of the effect of experimental cuts on the pseudorapidity and 
transverse momentum of the final state charged lepton and missing energy of the neutrino. 

The high luminosity (10 33 — 10 34 cm 2 s _1 ) at the LHC insures that systematic errors 
will play a dominant role in determining the accuracy of the cross-section. Thus, we present 
an analysis of the effect of the theoretical uncertainty in the evaluation of the acceptance, 
and propose which among the various available MC generators and computational schemes 
should be implemented to best simulate the events. 
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This paper is organized as follows. Sec. || will give an overview of the calculation and 
the computational tools used in the analysis. The next four sections are each devoted to 
estimating a class of systematic errors: electroweak corrections in Sec. ||, NNLO QCD in 
Sec. ||, QCD scale dependence in Sec. [|, and parton distribution function uncertainties in 
Sec. []. Finally, the results are compiled and summarized in Sec. |7| 

2. Theoretical Calculations and MC Generators 

The dominant production mechanism for Z or W bosons is the Drell-Yan process ]|], in 
which a quark and antiquark annihilate to form a vector boson, which subsequently decays 
into a lepton pair. The W production process is actually observed through the charged 
lepton and missing energy of the neutrino produced in its decay. In general, the cross- 
section may be inferred from the number Nyy B of observed events via the relation 

N$ s = a tot BR(W -» lv)A w j Cdt. (2.1) 

Ay/ is the acceptance obtained after applying the experimental selection criteria. For 
example, if the cuts require p T > p™ m , < rf < rj max , and^ T >^™ in , then 

1 rV~s/2 ^ /-Vi/2 

O-tot BR(W -> iv) J p ^in d?T Jymin 

rVmax r°° d 4 a 

/ dm dT]v dr> * dv "dvedr, BR( - W ~+ lv) ' (2 ' 2) 

Alternatively to the W production cross-section measurement, the corrected W yield 
can be used as a standard candle for a luminosity monitor in LHC if one calculates the 
cross-section and solves for J Cdt. The theoretical cross-section may be constructed by 
convoluting a parton-level cross-section for partons a and b with the parton density 
functions (PDFs) f a , fb for these partons, 

a th BR(W -^ev)=J2 f ^l^2/a(xi)/ 6 (x 2 ) a ab (xi,x 2 ), (2.3) 
a,b Jo 

integrating over the momentum fractions x\,X2, and applying cuts relevant to the exper- 
iment. Theoretical errors come from limitations in the order of the calculation of o a b, on 
its completeness (for example, on whether it includes electroweak corrections, on whether 
any phase space variables or spins have been averaged), and from errors in the PDFs. 

Since the final state may include additional partons which form a shower, the output 
from the hard QCD process must be fed to a shower generator to generate a realistic final 
state seen in a detector. This is possible only if the cross-section is simulated in an event 
generator. Calculating the acceptance for all but the simplest cuts will normally require 
an event generator as well. 

Thus, when constructing a simulation of an experiment, there is a range of choices 
which can be made among the tools currently available. An efficient calculation requires 
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selecting those adequate to meet the anticipated precision requirements, without perform- 
ing unnecessarily complex calculations. For example, while NNLO calculations are now 
available, the cross-sections are very complicated, do not always converge well, and require 
substantial time to calculate. For certain choices of cuts, it may be found that the effect of 
the NNLO result can be minimized, or that it can be represented by a simplified function 
for the parameters of interest. We will compare several possible schemes for calculating 
the W production cross-section and acceptance, and consider the systematic errors arising 
for these schemes. 

The most basic way to generate events is through one of the showering programs, such 
as PYTHIA [|, HERWIG @, IS A JET |l| or SHERPA @. These vary somewhat in 
their assumptions and range of effects included, but they all start with hard partons at a 
high energy scale and branch to form partons at lower scales, which permits a description 
of hadronization and realistic events. On their own, these programs typically rely on a 
leading order hard matrix element and include only a leading-log resummation of soft 
and collinear radiation in the shower, limiting their value in describing events with large 
transverse momentum. In addition, ISAJET lacks color-coherence, which is important in 



predicting the correct distribution of soft jets [13]. 

Fully exclusive NLO QCD calculations are available for W and Z boson production |U 
The MC generator MC@NLO [jl5| combines a parton-level NLO QCD calculation with the 
HERWIG ill parton shower, thus removing some of the limitations of a showering program 
alone. 

Since a sa a 2 s at LHC energies, NLO electroweak (EWK) corrections should appear at 
the same order as NNLO QCD. The MC@NLO package is missing EWK corrections, but 
the contribution of final state radiation (FSR) can be obtained by combining MC@NLO 
with PHOTOS [fig ], an add-on program which generates multi-photon emission from events 
created by the host program. In the case of Z exchange, FSR was expected to be the 
dominant contribution. [17] Another program, HORACE pj| , is available which includes 



exact 0{a) EWK corrections together with a final state QED parton shower. It was 
confirmed ||] that PHOTOS and HORACE agree within 1% for Z production. However, 
for W production, the W itself can radiate, so we should expect that the agreement may 
not be so close in this case. Recent studies [jO]] have emphasized the importance of the 
interplay between QCD and EWK corrections, particularly in the high p T tail important 
for new physics searches. 

The other available NLO and NNLO calculations are implemented as MC integra- 
tions, which can calculate a cross-section but do not provide unweighted events. Some of 
these are more differential than others. For example, the NNLO rapidity distribution is 
available in a program Vrap ||, but this distribution alone is not sufficient to calculate 
acceptances with cuts on the lepton pseudorapidities and transverse momenta. A differ- 
ential version of this NNLO calculation is implemented in a program FEWZ §, but this 
is not an event generator. Another available program is ResBos-A pip| , which resums soft 
and collinear initial state QCD radiation to all orders and includes NLO final state QED 
radiation. Resummation is expected to have advantages in realistically describing the small 
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Our analysis is conducted for di-lepton final states. The available calculations typically 
set the lepton masses to zero, so the lepton masses will be neglected throughout this paper 
and the choice of final state lepton has no effect on the calculations. In all results, t may be 
interpreted as either an electron or muon, v the appropriate accompanying (anti)neutrino, 
and we consider the case of production separately. We have chosen three sets of 
experimental cuts to reflect detector capabilities and to demonstrate the impact of physics 
effects on the acceptances depending on the selection criteria. 

3. Electroweak Corrections 

As noted above, both NNLO QCD corrections and NLO electroweak (EWK) corrections 
are expected to be needed to reach precisions on the order of 1% or better in W boson 



production. NLO electroweak^] and QCD corrections [14] are known both for W boson 
production and Z boson production. However, current state-of-the-art MC generators 
do not include both sets of corrections. The generator MC@NLO JO]] combines a MC 
event generator with NLO calculations of rates for QCD processes and uses the HERWIG 
event generator for the parton showering, but it does not include EWK corrections. Final 
state QED can be added using PHOTOS |l6jl , a process-independent module for adding 
multi-photon emission to events created by a host generator. However, some 0(a) EWK 
corrections are still missing, in particular radiation from the charged W itself, which was 
not an issue in Z production. 



To study the error arising from missing 0(a) EWK corrections, we used HORACE [18], 
a MC event generator that includes initial and final-state QED radiation in a photon shower 
approximation and exact 0{a) EWK corrections matched to a leading-log QED shower. 
To determine the magnitude of the error, we then compared the results from this generator 
to a Born- level calculation with final-state QED corrections added by PHOTOS. 

Specifically, we compared pp — > W — > Iv events generated by HORACE with the full 
0(a) corrections and parton-showered with HERWIG, to these events generated again by 
HORACE, but without EWK corrections (Born-level), showered with HERWIG+PHOTOS. 
CTEQ6.5M parton distribution functions p2| were used in the calculations. The results 
are shown in Tables ^ - || and in Figs. [I] - 1|[ 

For Tables || - ||, a standard set of cuts is used. We choose three sets of experimental 
cuts described in Table || to reflect detector capabilities and to demonstrate the impact of 
physics effects on the acceptances depending on the selection criteria. Here, r\ and p T are 
the pseudorapidity and transverse momentum of the final state charged leptons, and^ T 
is the missing transverse energy carried by the neutrino. The different cuts provide useful 
separation for between regions of the W spectrum which have different sensitivities to some 
of the sources of uncertainties. 

The first row in the table shows the total generator-level cross-sections before QCD parton 
showering (identified by the label "Total"). The Born+FSR column shows the effect of 
applying final state radiation (FSR) corrections only via PHOTOS. In PHOTOS, FSR 
affects the rates through the cuts only. 
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Transverse 




Missing Transverse 




Momentum (GeV/c) 


Pseudorapidity 


Energy (GeV) 


Cut 1 


p r > 25 


\v\ < i 


# T > 20 


Cut 2 


p T > 25 


K H < 2.2 


# T > 20 


Cut 3 


p r > 25 


H<i 


# T > 30 



Table 1: Acceptance regions for the Electroweak and NNLO studies 



Photonic and Electroweak Corrections: W + Production 





Born 


Born±FSR 


ElectroWeak 


Difference 


a (Total) 
a (Cut 1) 
a (Cut 2) 
a (Cut 3) 
A (Cut 1) 
A (Cut 2) 
A (Cut 3) 


10780.6 ± 1.2 
1321.3 ±5.4 
1639.1 ±5.9 
1008.9 ±4.9 
0.1652 ±0.0007 
0.2049 ± 0.0007 
0.1261 ±0.0006 


10780.6 ± 1.2 
1294.0 ± 5.4 
1604.7 ±5.9 
989.0 ± 4.8 
0.1618 ±0.0007 
0.2006 ± 0.0007 
0.1236 ±0.0006 


11201.4 ±1.6 
1345.8 ± 5.6 
1656.7 ±6.1 
1023.1 ±5.0 
0.1619 ±0.0007 
0.1993 ±0.0007 
0.1231 ±0.0006 


3.90 ± 0.02% 
4.00 ± 0.61% 
3.24 ± 0.53% 
3.45 ± 0.71% 
0.09 ± 0.59% 
-0.64 ±0.51% 
-0.44 ± 0.69% 



Table 2: Calculation of the W + — > £ + V£ cross-section a and acceptance A for various EWK 
corrections generated using HORACE 3.1, for £ = e or fi. 



Photonic and Electroweak Corrections: W Production 





Born 


Born±FSR 


ElectroWeak 


Difference 


a (Total) 
a (Cut 1) 
a (Cut 2) 
a (Cut 3) 
A (Cut 1) 
A (Cut 2) 
A (Cut 3) 


7998.1 ±1.1 
1498.5 ± 5.7 
1611.0 ±5.9 
1158.9 ±5.1 
0.1873 ±0.0007 
0.2014 ±0.0007 
0.1449 ±0.0006 


7998.1 ±1.1 

1466.7 ±5.7 

1578.8 ±5.8 
1135.4 ±5.1 

0.1833 ±0.0007 
0.1974 ±0.0007 
0.1419 ± 0.0006 


8311.3 ±2.6 
1524.1 ±5.9 

1648.5 ±6.1 

1180.6 ±5.3 
0.1834 ±0.0007 
0.1983 ±0.0007 
0.1420 ±0.0006 


3.92 ± 0.04% 
3.91 ± 0.57% 
4.41 ± 0.54% 
3.98 ± 0.66% 
0.05 ± 0.54% 
0.46 ± 0.50% 
0.07 ± 0.60% 



Table 3: Calculation of the W — > I vi cross-section a and acceptance A for various EWK 
corrections generated using HORACE 3.1, for £ = e or fi. 

The ElectroWeak column includes the full HORACE EWK corrections. In the final 
column, we give the difference between the previous two columns, to compare the full 
EWK correction FSR alone. The cross-section (acceptance) results show agreement within 
4% (0.7 %) between the two schemes. The maximum error in the cross-section is 4.4% 
corresponding to cut 2 for W~ or 4% corresponding to cut 1 for W + , and the maximum 
error in the acceptance is 0.64 % or 0.46 % for the cut 2 for W + or W~ correspondingly. 

4. NNLO QCD Uncertainties 

QCD uncertainties include errors due to missing higher-order corrections in the hard matrix 
element, uncertainties in the parton distribution functions, and approximations made in the 
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Generated W* -> Tv/ny) M x Generated W + -> Tv, M,. 




20 40 60 80 100 120 140 160 20 40 60 80 100 120 140 160 



l^ w ,(GeV/c 2 ) l^ w ,(GeV/c 2 ) 

(a) (b) 

Figure 1: Comparison of the generated (a) W + boson transverse mass distributions and (b) i + v^ 
transverse mass distributions for the process W + — > t + vi(n^) in HORACE 3.1 including 0(a) 
EWK corrections showered with HERWIG (open red squares), HORACE Born-level showered with 
HERWIG plus PHOTOS (black circles), and HORACE Born-level (blue stars). 



Generated W" -» I Wrryj M x Generated W" -> IV M T 




20 40 60 80 100 120 140 160 20 40 60 80 100 120 140 160 



M TW (GeV/c 2 ) M TW (GeV/c 2 ) 

(a) (b) 

Figure 2: Comparison of the generated (a) W~ boson transverse mass distributions and (b) l~Vi 
transverse mass distributions for the process W~ — > l~U^(rvy) in HORACE 3.1 including O(a) 
EWK corrections showered with HERWIG (open red squares), HORACE Born-level showered with 
HERWIG plus PHOTOS (black circles), and HORACE Born-level (blue stars). 



showering algorithms. In the following, we will evaluate the errors introduced by omitting 
the NNLO corrections by using MC@NLO, and calculate K- factors which can be used to 
introduce NNLO corrections to the MC@NLO calculation. We will also examine the effect 
of uncertainties in the PDFs. For these studies, we choose the same three sets of cuts as 
given previously in Table |]. 

We begin by examining the NNLO corrections, using the state-of-the-art program 
FEWZ Q , which is differential in the and the lepton transverse momenta and pseudorapidi- 
ties. The FEWZ program is at NNLO in perturbative QCD, includes spin correlations, and 
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Number of FSR Photons 



FSR Photon Transverse Momentum 



HOR Born+Herwig 

HOR Born+Herwig+PHOTOS 

HOR EWK+Herwig 








z "5" 




eI 








E« 10 




z 




z a. 




z ^- 




E oT 




- "O 













HOR Born+Herwig+PHOTOS 
HOR EWK+Herwig 




Figure 3: Comparison of the number n of 
hard final state radiation (FSR) photons in 
W -> lv{ni) for HORACE 3.1 including 
0(a) EWK corrections showered with HER- 
WIG (open red squares), HORACE Born- 
level showered with HERWIG plus PHO- 
TOS (black circles), and HORACE Born- 
level (blue stars). 



2 4 6 8 10 12 14 

P tY (GeV/c) 

Figure 4: Comparison of W — ► lv(n~j) fi- 
nal state radiation (FSR) transverse momen- 
tum distributions for HORACE 3.1 including 
EWK corrections showered with HERWIG 
(open red squares) and HORACE Born- 
level showered with HERWIG plus PHOTOS 
(black circles). 



takes into account finite widths effects. Since we are interested primarily in studies about 
the W peak, we choose the renormalization and factorization scales to be n F = fi R = Mw- 
Scale dependence will be discussed in detail the next section. 

A comparison of the effects of higher order QCD corrections on the cross-section and 
acceptance is presented in Tables f| - |5] and Figs. || - ||. Both the NLO and NNLO calcula- 
tions are done with CTEQ6.5M PDFs |22|], since we will be calculating X-factors intended 
to rescale NLO calculations which have used these PDFs. However, since CTEQ6.5M 
PDFs are only available at NLO, we have also repeated the calculations in the tables using 
MRST PDFs, and found the results to be compatible, as will be discussed later in this 



NNLO Cross Sections a (pb) for W + Production 



Cut 


MC@NLO 


FEWZ NLO 


FEWZ NNLO 


K- fact or 


Kx MC@NLO 


a tot 
1 

2 
3 


12587 ± 13 
2072.6 ±6.0 
2541.8 ±6.5 
1575.0 ± 5.4 


12869 ± 12 
2157.6 ±2.2 
2682.4 ± 2.7 
1656.0 ± 1.6 


12780 ± 48 
2122.9 ± 14.7 
2651.6 ±21.2 
1605.8 ± 13.1 


0.9931 ± 0.0039 
0.9839 ± 0.0069 
0.9885 ± 0.0079 
0.9970 ± 0.0080 


12500 ± 51 

2039.2 ± 15.5 
2512.6 ±21.1 

1570.3 ± 13.7 


NNLO Acceptances (%) for W + Production 


Cut 


MC@NLO 


FEWZ NLO 


FEWZ NNLO 


K- fact or 


Kx MC@NLO 


1 

2 
3 


16.47 ±0.05 
20.19 ±0.05 
12.51 ±0.04 


16.77 ±0.02 
20.84 ±0.03 
12.87 ±0.02 


16.61 ±0.13 
20.75 ±0.18 
12.57 ±0.11 


0.9908 ± 0.0079 
0.9954 ± 0.0089 
0.9765 ± 0.0089 


16.31 ±0.14 
20.10 ±0.19 
12.22 ±0.12 



Table 4: Calculation of the W + — > i+vt (£ = e or fx) cross-section at NLO using MC@NLO, and 
at NLO and NNLO using FEWZ, for the cut region defined in Table |[ 
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Lepton p T Cut: (M > 20, \t\\ < 1 .0) 




15 20 25 30 35 
P Tmi „( GeV/C > 

(a) 



Lepton p T Cut: (M > 20, ]r\\ < 1 .0) 




P T ■ (GeV/c) 

"T mm v ' 



(b) 



Figure 5: Cross-section (a) and acceptance (b) versus cut on lepton p T at NLO (hashed bands) 
and NNLO (points), as calculated using FEWZ for W + -> l + v t . 



section. All results in the tables are calculated at scale Myy. In the figures, the NLO 
results are displayed as a band spanning the range of scales from to 2M\y. The 

scale dependence of the NNLO result is small enough to be comparable to the precision of 
the MC evaluation of the integrals, so only the average of the high and low scales is plotted, 
with error bars reflecting a combination of statistical and scale variation uncertainties. 

Since the NNLO matrix element has not yet been interfaced to a shower, we cannot 
directly compare FEWZ to MC@NLO. The best we can do at this time is to use MC@NLO 
to obtain the NLO showered result, and multiply this by a iT-factor obtained by taking 
the ratio of the NNLO to NLO results derived from FEWZ. This procedure is reasonable 
except in threshold regimes where the fixed-order NLO result in FEWZ would be unreliable. 
This is similar to methods that have been used for calculating NNLO corrections to Higgs 
production |2^]. The differences of these K- factors from unity are shown in Figs. |9| and 10 
for both the cross-sections and acceptances, as a function of cuts on the lepton p T and r\. 



NNLO Cross Sections a (pb) for W Production 



Cut 


MC@NLO 


FEWZ NLO 


FEWZ NNLO 


K- fact or 


Kx MC@NLO 


a tot 
1 

2 
3 


9202.5 ±9.9 
1733.5 ±4.6 
1858.3 ±4.8 
1341.2 ±4.2 


9450.5 ±9.2 

1794.0 ± 1.8 
1950.7 ±1.9 

1404.1 ± 1.4 


9357.6 ± 34.7 
1772.1 ± 14.3 
1890.4 ±21.4 
1355.1 ± 13.1 


0.9902 ± 0.0038 
0.9878 ± 0.0080 
0.9691 ±0.0110 
0.9651 ± 0.0094 


9112.3 ±36.3 

1712.4 ± 14.6 
1800.9 ±21.0 
1294.4 ± 13.2 


NNLO Acceptances (%) for W Production 


Cut 


MC@NLO 


FEWZ NLO 


FEWZ NNLO 


K- fact or 


Kx MC@NLO 


1 

2 
3 


18.84 ± 0.05 
20.19 ±0.05 
14.58 ± 0.05 


18.98 ± 0.03 
20.64 ±0.03 
14.86 ± 0.02 


18.94 ±0.17 
20.20 ± 0.24 
14.48 ±0.15 


0.9976 ± 0.0090 
0.9787 ±0.0118 
0.9747 ± 0.0102 


18.79 ±0.18 
19.76 ± 0.24 
14.21 ±0.16 



Table 5: Calculation of the W — > £ vg, (£ = e or fx) cross-section at NLO using MC@NLO, and 
at NLO and NNLO using FEWZ, for the cut region defined in Table |[ 
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Lepton |ti| Cut: (M > 20, p y > 25) 



Lepton (pi Cut: (M > 20, p T > 25) 




(a) 




I Lip- 



(b) 



Figure 6: Cross-section (a) and acceptance (b) versus cut on lepton |ry| at NLO (hashed bands) 
and NNLO (points), as calculated using FEWZ for W + -»■ l + v t . 



Lepton p T Cut: (M > 20, (nl < 1 .0) 



Lepton p Cut: (M > 20, |r|| < 1 .0) 
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p T . (GeV/c) 
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(b) 



Figure 7: Cross-section (a) and acceptance (b) versus cut on lepton p T at NLO (hashed bands) 
and NNLO (points), as calculated using FEWZ for W~ — ► l~vi. 



The resulting accepted cross-section is shown in the Kx MC@NLO column of Tables Q - 
|5] and in Figs. 11 and 12 as a function of the same cuts as in Figs. ^ and |l^. The size of 
K — 1 is a good indicator of the error due to missing NNLO if MC@NLO is used without 
corrections. 

The results in Tables |] - || show i^T-factors corresponding to an NNLO correction of 
about 1% for the cross-sections, or up to 3.5% for the Cut 3 cross-section. NNLO 
corrections to the acceptances were generally at the 2.5% level or less, with Cut 3 again 
having the biggest corrections. Cut 1 acceptances had a 1% or less correction for both W^. 
Convergence of he FEWZ W calculation appears better than for the Z calculation, giving 
more rapid convergence and smaller Monte-Carlo errors in comparable run-times. Thus, 
1% evaluations were generally possible in the W calculation, but were not always practical 
in the corresponding Z calculations. || 

As noted above, the CTEQ6.5M PDFs do not include NNLO corrections, it is impor- 
tant to check their compatibility with the NNLO calculation by repeating the calculation 
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Figure 8: Cross-section (a) and acceptance (b) versus cut on lepton \rj\ at NLO (hashed bands) 
and NNLO (points), as calculated using FEWZ for W~ — > l~T>i. 



if -factors with CTEQ and MRST PDFs 



Cut 


CTEQ K(a) 


MRST K(a) 


CTEQ K(A) 


MRST K(A) 


W + cut 1 
W~ cut 1 
W + cut 2 
W~ cut 2 
W+ cut 3 
W~ cut 3 


0.9839 ± 0.0069 
0.9878 ± 0.0080 
0.9885 ± 0.0079 
0.9691 ±0.0110 
0.9697 ± 0.0080 
0.9651 ± 0.0094 


0.9767 ± 0.0081 
0.9538 ± 0.0086 
0.9745 ± 0.0100 
0.9712 ±0.0109 
0.9510 ±0.0097 
0.9430 ± 0.0209 


0.9908 ± 0.0079 
0.9976 ± 0.0090 
0.9954 ± 0.0089 
0.9787 ±0.0118 
0.9765 ± 0.0089 
0.9747 ±0.0102 


0.9988 ± 0.0091 
0.9804 ± 0.0097 
0.9966 ±0.0110 
0.9984 ±0.0119 
0.9725 ±0.0106 
0.9694 ±0.0218 



Table 6: Comparison of the if factors for cross-sections and acceptances calculated using 
CTEQ6.5M (NLO) and MRST 2002 (NNLO) PDFs for both W+ and W~ production. K{a) 
is the if- factor for the cross-section and if (A) is the if- factor for the acceptance, as displayed in 
Tables | - 1. 



with MRST 2002 NNLO PDFs. It is not obvious that substituting NLO PDFs in this cal- 
culation is reasonable, but in fact, it turns out to be. The if -factors for the cross-section 
and acceptance for each of the three cuts is shown in Table || for both sets of PDFs. The 
if (<r)-factors are ratios of NNLO to NLO cross-sections, and if {A) are ratios of NNLO 
to NLO acceptances relative to the total cross section. The agreement of the if factors 
for the cross sections varies from less than 1% to 3.6%, with the maximum disagreement 
for W~ cut 1. The acceptances agree to within 1% except for W~ cuts 1 and 2, where 
the agreement is within 2%. There is a still some variation in the results due to the slow 
convergence of the NNLO calculation, typically in the 1 — 2% range. Only cut 1 for the 
W~ has a PDF dependence outside that range, and this could still be primarily due to 
Monte Carlo errors. 



5. Scale Dependence 

Perturbative QCD calculations at fixed order depend on the factorization and renormal- 
ization scales introduced in the calculation. Thus, the previous calculations have an added 
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Figure 9: Fractional difference in the NNLO and NLO cross-sections (left-hand side) and accep- 
tances (right-hand side) as a function of the lepton (a) p , and (b) \r)\ cuts as in Figs. Figs || - ^ 
for W + — > l + vt. These differences are the factor K — 1 for the cross-section and acceptances, 
respectively. 



uncertainty due to the choice of certain fictitious scales appearing in the calculation. In 
a complete, all order calculation, or one completely resummed in the soft and collinear 
regimes and properly matched to the PDFs, there would be no dependence on these scales. 
However, in a fixed order calculation matched to PDFs, a dependence on the factorization 
scale fjL F and renormalization scale fi R appear in the final results. The effect of scale choice 
is significant at NLO. Adding NNLO effects is found to reduce scale dependence consider- 
ably |,[7|, though it can remain significant near thresholds where NLO effectively becomes 
leading order. 

As is customary, we will choose the renormalization and factorization scales to be 
identical, and investigate the scale dependence by varying them by a factor of 2 or 1/2 
about a central value of n F R = Mz, which is typical of the scales in our acceptance, and 
was the central value chosen in the previous section. 

Tables |I| - |B| included only the central scale My/ . Tables show the total cross 

sections and acceptances for lepton production calculated by FEWZ at three different 
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Figure 10: Fractional difference in the NNLO and NLO cross-sections (left-hand side) and ac- 
ceptances (right-hand side) as a function of the lepton (a) p and (b) |rj| cuts as in Figs. Figs fj] 
- H for W~ — » l~vt. These differences are the factor K — 1 for the cross-section and acceptances, 
respectively. 



renormalization and factorization scales Mw/2, Myy, and 2M\y. The acceptances for 
the final state leptons are as defined in Table [l]. For a measure of the size of the scale 
dependence, the final column of each table shows the maximum difference between the 
three values divided by average, with an error calculated assuming the statistical errors in 
the three MC runs in each row are independent. 

We can see that the scale dependence of the cross-sections at NLO is typically of 
order ±6%. The scale dependence of NLO acceptances is dramatically reduced due to 
correlations in the scale dependence of the cut and uncut cross-sections used to compute 
it. Adding NNLO reduces the scale dependence of the cross-sections to less than 2% in 
most cases. For some cuts, the convergence of the Monte-Carlo is a significant limitation 
on the accuracy of this result, since the NNLO MC errors can be comparable to the scale 
dependence. 
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Figure 11: Accepted NNLO cross-section estimated from MCQNLO scaled by the if -factor versus 
the lepton (a) p T , and (b) cuts for W + — > t + V£. 
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Figure 12: Accepted NNLO cross-section estimated from MC@NLO scaled by the if -factor versus 
the lepton (a) p and (b) \r)\ cuts for W~ — > £~ De- 



fy. Uncertainties Due to the Parton Distribution Function 

Phenomenological parameterizations of the PDFs are taken from a global fit to data. There- 
fore, uncertainties on the PDFs arising from diverse experimental and theoretical sources 
will propagate from the global analysis into the predictions for the W/Z cross-sections. 



Figure 12 and Figure 14 show the results of the inclusive W to di-lepton production cross- 
section using various CTEQ |§ and MRST (2§ PDFs. The upward shift of about 7% 
(between CTEQ6.1 and 6.5 and MRST2004 and 2006) results from the inclusion of heavy 
quark effects in the latest PDF calculations. The acceptance due to the cuts in Table [j] 
using each of these PDFs is shown in Fig. ^ and Fig. 16. 1 

The uncertainties in the PDFs arising from the experimental statistical and systematic 
uncertainties, and the effect on the production cross-section of the W boson, have been 



1 Theoretical issues which may affect the contribution of the PDFs to the NNLO if- factor are not included 
here as we are concerned primarily with the error at NLO. See Refs. fed] for details. 
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Total Cross-Section (in pb): W + 



Order 


M w /2 M w 2M W 


Aa/a 


NLO 
NNLO 


12517 ± 11.6 12869 ± 11.9 13203 ± 13.2 
12708 ± 63.6 12780 ± 48.2 12903 ± 49.7 


0.0534 ± 0.0013 
0.0153 ± 0.0060 


Cut Region 1: W+ 


Order 


M w /2 M w 2M W 


Act /a 


NLO 
NNLO 


2096.9 ±2.1 2157.6 ± 2.2 2220.1 ±2.2 
2123.3 ±23.2 2122.9 ± 14.7 2139.8 ± 15.7 


0.0571 ±0.0014 
0.0079 ± 0.0121 


Cut Region 2: W + 


Order 


M w /2 M w 2M W 


Aa/a 


NLO 
NNLO 


2608.3 ± 2.6 2682.4 ± 2.7 2762.0 ± 2.7 
2606.5 ± 24.6 2651.6 ±21.2 2643.6 ± 19.5 


0.0572 ± 0.0014 
0.0172 ±0.0117 


Cut Region 3: W + 


Order 


M w /2 M w 2M W 


Aa/a 


NLO 
NNLO 


1607.6 ± 1.6 1656.0 ± 1.6 1703.3 ± 1.7 
1620.4 ± 11.8 1605.8 ± 13.1 1625.6 ± 13.7 


0.0578 ± 0.0014 
0.0122 ±0.0113 



Table 7: Scale dependence of the total and accepted cross-sections (in pb) for W + boson production 
calculated by the FEWZ program at order NLO and NNLO. The final column is a measure of scale 
dependence obtained by dividing the maximum spread by the average for the three points. 



Cut Region 1: W + (% Accepted) 



Order 


M w /2 


M w 


2M W 


AA/A 


NLO 


16.75 ±0.02 


16.77 ±0.02 


16.82 ±0.02 


0.0037 ± 0.0020 


NNLO 


16.71 ±0.20 


16.61 ±0.13 


16.58 ±0.14 


0.0076 ±0.0136 




Cut Region 2: W + 


(% Accepted) 




Order 


M w /2 


M w 


2M W 


AA/A 


NLO 


20.84 ±0.03 


20.84 ±0.03 


20.92 ±0.03 


0.0039 ± 0.0020 


NNLO 


20.51 ±0.22 


20.75 ±0.18 


20.49 ±0.17 


0.0127 ±0.0132 




Cut Region 3: W + 


(% Accepted) 




Order 


M w /2 


M w 


2M W 


AA/A 


NLO 


12.84 ±0.02 


12.87 ±0.02 


12.90 ±0.02 


0.0045 ± 0.0020 


NNLO 


12.75 ±0.11 


12.57 ±0.11 


12.60 ±0.12 


0.0147 ±0.0128 



Table 8: Scale dependence of the acceptances A in the various cut regions for W + boson production 
calculated by the FEWZ program at order NLO and NNLO. The final column is a measure of scale 
dependence obtained by dividing the maximum spread of the three preceding columns by their 
average. 



studied using the standard methods proposed in Refs. [22, 25 1. For the standard set of 
PDFs, corresponding to the minimum in the PDF parameter space, a complete set of 
eigenvector PDF sets, which characterize the region nearby the minimum and quantify its 
error, have been simultaneously calculated. From the minimum set and these "error" sets 
we calculate the best estimate and the uncertainty for the W cross-section. We do this 
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Total Cross-Section (in pb): W 



Order 


M w /2 M w 2M W 


Aa/a 


NLO 
NNLO 


9184.6 ±9.1 9450.5 ±9.2 9706.5 ±9.5 
9307.3 ± 34.3 9357.6 ± 34.7 9419.5 ± 36.0 


0.0553 ± 0.0014 
0.0120 ± 0.0053 


Cut Region 1: W~ 


Order 


M w /2 M w 2M W 


Act jo 


NLO 
NNLO 


1741.4 ± 1.7 1794.0 ± 1.8 1842.4 ± 1.8 

1717.5 ± 13.5 1772.1 ± 14.3 1767.8 ± 12.6 


0.0563 ± 0.0014 
0.0312 ± 0.0109 


Cut Region 2: W~ 


Order 


M w /2 M w 2M W 


Act J a 


NLO 
NNLO 


1894.0 ± 1.9 1950.7 ± 1.9 2005.5 ± 2.0 
1892.9 ± 15.7 1890.4 ±21.4 1906.4 ± 15.0 


0.0572 ± 0.0014 
0.0084 ± 0.0131 


Cut Region 3: W~ 


Order 


M w /2 M w 2M W 


Act /a 


NLO 
NNLO 


1360.2 ± 1.4 1404.1 ± 1.4 1449.3 ± 1.5 
1343.8 ± 11.8 1355.1 ± 13.1 1358.7 ±13.7 


0.0634 ± 0.0014 
0.0110 ±0.0135 



Table 9: Scale dependence of the total and accepted cross-sections (in pb) for W~ boson production 
calculated by the FEWZ program at order NLO and NNLO. The final column is a measure of scale 
dependence obtained by dividing the maximum spread by the average for the three points. 



Cut Region 1: W~ (% Accepted) 



Order 


M w /2 


M w 


2M W 


AA/A 


NLO 


18.96 ± 0.03 


18.98 ±0.03 


18.98 ±0.03 


0.0013 ± 0.0020 


NNLO 


18.45 ±0.16 


18.94 ±0.17 


18.77 ±0.15 


0.0259 ±0.0121 




Cut Region 2: W + 


(% Accepted) 




Order 


M w /2 


M w 


2M W 


A A/ A 


NLO 


20.62 ± 0.03 


20.64 ±0.03 


20.66 ± 0.03 


0.0019 ±0.0020 


NNLO 


20.34 ±0.18 


20.20 ±0.24 


20.24 ±0.18 


0.0067 ± 0.0142 




Cut Region 3: W + 


(% Accepted) 




Order 


M w /2 


M w 


2M W 


AA/A 


NLO 


14.81 ±0.02 


14.86 ±0.02 


14.93 ±0.02 


0.0082 ± 0.0020 


NNLO 


14.44 ±0.14 


14.48 ±0.15 


14.42 ±0.16 


0.0039 ± 0.0145 



Table 10: Scale dependence of the acceptances A in the various cut regions for W~ boson produc- 
tion calculated by the FEWZ program at order NLO and NNLO. The final column is a measure 
of scale dependence obtained by dividing the maximum spread of the three preceding columns by 
their average. 



using the Hessian error method, where the cross-section results from the various eigenvector 
PDF sets have been combined according to the prescriptions found in [p2j ]. Fig. 13 and 
Fig. |i~4}list the results for the different PDFs and Table 11 summarizes the results of the 



latest CTEQ and MRST PDF sets. The difference in the uncertainties (approximately a 
factor of two) between the results obtained from the CTEQ and MRST PDF error sets is 
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Figure 13: Comparison of W —*lu total cross-sections for several recent PDF calculations. 
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Figure 14: Comparison of W + — > l + v total cross-sections for several recent PDF calculations. 



due to different assumptions made by the groups while creating the eigenvector PDF sets. 
Finally we study the sensitivity of the kinematic acceptance calculations to the un- 



certainties affecting the PDF sets. Figs. O [19| for W and Figs. [18] EUI for W + show the 



systematic error on the production cross-sections as a function of the \rj\ cut and minimum 
lepton p T for variations on the three types of cuts in Table |l[ The fractional uncertainties, 
shown in in the same figures, demonstrate that the relative uncertainty in the cross-section 
is very flat as a function of the kinematic cuts, until the region of extreme cuts and low 
statistics in the MC are reached. The corresponding uncertainty on the acceptance as a 
function of the kinematic cuts is shown in Figs ^ and 23 for W~ and Figs 22 and 24 
for W + . These show a similar dependence to the cross-section uncertainties, though the 
fractional errors are smaller. 



7. Conclusions 



To evaluate the overall contribution from theoretical uncertainties to both the cross-section 
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and acceptance calculations for the decay mode W — > t vg (£ = e or //) at the LHC we add 
the uncertainties from each of the sources considered in the preceding sections. We compile 
the errors assuming that the calculation is done with MC@NLO at scale fi F = fi R = M\y 
and interfaced to PHOTOS to add final state QED radiation. The missing electroweak 
contribution may then be inferred from HORACE as in Sec. |3|. For these errors we take 
those resulting from the tight cut set in Tables Q - |3| The cuts described in Table |l| are 
considered the most representative of likely analysis cuts for the LHC experiments. 

QCD uncertainties may be divided into two main classes. If the NNLO if-factor is set 
to 1, there is a missing NNLO contribution <5nnlo = K — 1. Since K has residual NNLO 
scale dependence, we must also take this into account and write K = 1 + ^nnlo ± ^ scale- 
The factor 5nnlo can be inferred from Tables Q - [5] and 5 scale can be inferred from half the 
scale variation of the NNLO entries in Tables 0-|9|, as discussed in detail in the conclusions 
of Ref. §. 

Both classes of QCD errors are also associated with a "technical precision" due to 
limitations of the computing tools used to evaluate them. Significant improvements in the 
NNLO precision could be obtained if a program with faster convergence were available. We 
therefore include an "error on the error" for the QCD errors, and propagate these through 
in the usual fashion to derive a final accuracy for the total QCD uncertainty estimate. 
This sets a limitation on how much the NLO calculation can realistically be improved 
using currently available NNLO results. 

These contributions to QCD errors are summarized in Table 12 for the total cross- 
section and the three cuts of Table |]. Results are shown both for the three cut cross-sections 
and their ratio for the total cross-section, and the errors are assumed to be uncorrelated. 

If the K- factor had not been calculated at NNLO, the error of the NLO cross-sections 
could have been roughly estimated from half the width of the scale-dependence band, or 
half the NLO results for Aa/a in Tables |?] and |9| , giving uncertainties of 2.8 — 3.2%. The 
errors calculated from the K factors are in within the limits these expectations, up to the 
technical precision of the calculation. A similar error NLO estimate for the error in the 
acceptance based on Tables |8| and |l(] would predict at most 0.4% missing NNLO. 

The final contribution to the total error considered here is the uncertainty from the 
PDFs. This may be extracted from the results of Sec. ^ by taking the errors from the 
CTEQ6.5 results for Cut 1 (see Table 0). The errors are asymmetric, so we take the largest 
of the two (up or down) uncertainties as the total fractional error for the PDF calculation. 
We choose the first cut set, since it is the most representative of likely analysis cuts at the 
LHC experiments. CTEQ errors, rather than the MRST errors, are used because they give 
a more conservative estimate. The difference between the results obtained by the latest 
CTEQ and MRST PDFs is less than the maximum error quoted for CTEQ for all three 
cut regions. 

The errors are added in quadrature, assuming no correlations, and the results are given 
in Table 13. There is, in fact, no concensus on the best way to combine these errors, so 
the total error should be considered an estimate. The QCD error is taken for Cut 1 for the 
same reasons as given above. In addition, as we have discussed above, we propagate the 
"error on the error" for each of the contributions in order to have some reasonable estimate 
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of the accuracy of the quoted total theoretical uncertainty. The exception to this is the 
PDF error, which can be considered as an upper limit on the uncertainty, and therefore 
does not need an additional accuracy. We conclude that the event generator MC@NLO 
interfaced to PHOTOS should be sufficient to guarantee an overall theoretical uncertainty 
on the W production due to higher order calculation, PDFs, and renormalization scale at 
the level of 5.5 — 5.9% for the total cross-section of W T and at approximately 2.5% for the 
acceptance of W^. 

W production will provide a valuable tool for studying QCD, measuring precision 
electroweak physics, and monitoring the luminosity. As the luminosity increases, the large 
statistics will permit a further improvement in the systematic uncertainties due to the 
PDFs. Adding complete 0(a) EWK corrections to the event generator would eliminate 
most of the EWK uncertainty, and incorporating NNLO QCD corrections would substan- 
tially reduce the QCD uncertainties. 

Reaching a combined precision of 1%, as desired in the later stages of analysis at high 
integrated luminosity, will require new tools. In addition to improved PDFs, an event gen- 
erator combining NNLO QCD with complete 0(a) EWK corrections will be needed, with 
exponentiation in appropriate regimes, and adequate convergence properties to technically 
reach the required precision. (2?]] Measurement strategies have been proposed that may mit- 
igate some of the effects of systematic errors on the precision of the W measurements.^] A 
combination of improved calculations and improved meaurements will be needed to permit 
the desired precision to be reached as the integrated luminosity increases to a point where 
it is needed. 
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Figure 15: Comparison of W~ — > l~V (£ - 
calculations for acceptance regions (a) Cut 1, 
left-hand plots show the total acceptance and 
acceptance. 



CTEQ6.5 MRST'06 (NNLO) GTEQ6.1 MRST'04 (NNLO)MRST04 (NLO) MRST'01 

PDFSet 

(c) 

= e or fj.) acceptances A, with several recent PDF 
(b) Cut 2, and (c) Cut 3, as defined in Table @. The 
the right hand plots show the fractional error on the 



- 19 - 





W~ 


W+ 




No cut 






Cut F 


Legion 1 






No cut 


Cut Region 1 


PDF Set 


a (pb) 


A<7 + 


Acr_ 


a (pb) 


Ao+ 


Acr_ 


A 


AA + 


AA_ 


a (pb) 


Aa + 


Acr_ 


a (pb) 


A<7 + 


Acr_ 


A 


AA+ 


AA_ 


CTEQ6.5 


9195 


501 


482 


1787.3 


57.5 


59.1 


0.185 


0.0028 


0.0041 


12560 


672 


691 


2107.26 


84.4 


83.5 


0.158 


0.0033 


0.0036 


MRST2006 


9206 


294 


344 


1824.2 


15.0 


67.0 


0.189 


0.0027 


0.0054 


12379 


459 


380 


2106.42 


23.7 


45.9 


0.160 


0.0010 


0.0031 


CTEQ6.1 


8552 


504 


503 


1671.1 


64.6 


98.9 


0.187 


0.0014 


0.0058 


11637 


763 


589 


1944.12 


115.6 


61.4 


0.157 


0.0047 


0.0028 


MRST2004 (NNLO) 


8702 


284 


292 


1681.8 


48.4 


22.7 


0.184 


0.0045 


0.0011 


11798 


397 


386 


1983.64 


56.9 


26.5 


0.159 


0.0028 


0.0024 


MRST2004 (NLO) 


8747 


286 


294 


1703.3 


49.1 


23.0 


0.186 


0.0045 


0.0011 


11655 


392 


382 


1949.68 


56.0 


26.1 


0.157 


0.0028 


0.0024 


1\ :fT> C T - -1(1(1 1 

MKo 1 2001 


8695 


284 


292 


1676.7 


48.3 


22.6 


0.184 


0.0045 


0.0011 


11743 


395 


385 


1937.29 


55.6 


25.9 


0.156 


0.0028 


0.0024 




No cut 


Cut Region 2 


No cut 


Cut Region 2 


PDF Set 


a (pb) 


A<7 + 


A<7_ 


a (pb) 


A<7 + 


A<7_ 


A 


AA+ 


AA_ 


a (pb) 


A(T + 


A(j_ 


a (pb) 


A(T + 


A(T_ 


A 


AA+ 


AA_ 


CTEQ6.5 


9195 


501 


482 


1907.6 


90.2 


64.5 


0.198 


0.0041 


0.0002 


12560 


672 


691 


2586.6 


81.1 


95.2 


0.193 


0.0018 


0.0026 


MRST2006 


9206 


294 


344 


1938.3 


38.1 


19.4 


0.201 


0.0035 


0.0001 


12379 


459 


380 


2568.4 


65.4 


13.8 


0.195 


0.0033 


0.0016 


CTEQ6.1 


8552 


504 


503 


1767.4 


95.6 


64.2 


0.197 


0.0049 


0.0012 


11637 


763 


589 


2379.47 


120.9 


75.9 


0.192 


0.0023 


0.0039 


n tt~» nrnaAr\ d /matt /~\ \ 

MRST2004 (NNLO) 


8702 


284 


292 


1835.2 


32.9 


31.6 


0.202 


0.0028 


0.0014 


11798 


397 


386 


2452.37 


49.3 


20.8 


0.195 


0.0019 


0.0015 


MRST2004 (NLO) 


8747 


286 


294 


1846.6 


33.1 


31.8 


0.201 


0.0028 


0.0014 


11655 


392 


382 


2420.16 


48.6 


20.5 


0.195 


0.0019 


0.0015 


MRST2001 


8695 


284 


292 


1831.9 


32.8 


31.5 


0.201 


0.0028 


0.0014 


11743 


395 


385 


2418.73 


48.6 


20.5 


0.194 


0.0019 


0.0015 




No cut 


Cut Region 3 


No cut 


Cut Region 3 


PDF Set 


a (pb) 


A(T + 


Acr_ 


a (pb) 


A(T + 


Acr_ 


A 


AA+ 


AA_ 


a (pb) 


A(T + 


Acr_ 


a (pb) 


A(T + 


Acr_ 


A 


AA+ 


AA- 


CTEQ6.5 


9195 


501 


482 


1362.5 


41.9 


51.8 


0.144 


0.0025 


0.0041 


12560 


672 


691 


1572.78 


73.0 


64.4 


0.120 


0.0031 


0.0030 


MRST2006 


9206 


294 


344 


1392.4 


15.8 


56.0 


0.147 


0.0006 


0.0047 


12379 


459 


380 


1572.68 


25.6 


37.8 


0.122 


0.0007 


0.0027 


CTEQ6.1 


8552 


504 


503 


1265.6 


65.5 


62.5 


0.144 


0.0019 


0.0034 


11637 


763 


589 


1451.63 


89.5 


43.9 


0.120 


0.0040 


0.0022 


MRST2004(NNLO) 


8702 


284 


292 


1279.6 


44.7 


15.1 


0.142 


0.0040 


0.0012 


11798 


397 


386 


1476.94 


32.5 


33.2 


0.121 


0.0015 


0.0031 


MRST2004 (NLO) 


8747 


286 


294 


1294.7 


45.3 


15.3 


0.144 


0.0041 


0.0012 


11655 


392 


382 


1451.63 


31.9 


32.6 


0.120 


0.0015 


0.0031 


MRST2001 


8695 


284 


292 


1273.1 


44.5 


15.0 


0.144 


0.0041 


0.0012 


11743 


395 


385 


1452.54 


32.0 


32.7 


0.119 


0.0015 


0.0031 



Table 11: Cross-sections a for , and acceptances A, with asymmetric Hessian uncertainties as calculated using several recent PDF sets for the 
three cut regions defined in Table 1. 



QCD Uncertainties (%) 



W + Cross-Section Act 


Uncertainty 


a tot 


Cut 1 


Cut 2 


Cut 3 


Missing NNLO 
Scale Dependence 


-0.69 ±0.39 
0.76 ±0.30 


-1.61 ±0.69 
0.40 ±0.61 


-1.15 ±0.79 
0.86 ±0.59 


-3.03 ± 0.80 
0.61 ±0.56 


Total 


1.03 ±0.34 


1.66 ±0.69 


1.43 ±0.73 


3.09 ±0.79 


Error in W + Acceptance (AA) 


Uncertainty 




Cut 1 


Cut 2 


Cut 3 


Missing NNLO 
Scale Dependence 




-0.92 ±0.79 
0.38 ± 0.68 


-0.46 ±0.89 
0.63 ±0.66 


-2.35 ±0.89 
0.74 ± 0.64 


Total 




1.00 ±0.78 


0.78 ± 0.75 


2.47 ±0.87 


W Cross-Section Ac- 


Uncertainty 


a tot 


Cut 1 


Cut 2 


Cut 3 


Missing NNLO 
Scale Dependence 


-0.98 ± 0.38 
0.60 ±0.26 


-1.22 ±0.80 
1.56 ±0.54 


-3.09 ± 1.10 
0.42 ± 0.66 


-3.49 ± 0.94 
0.55 ±0.68 


Total 


1.15 ±0.35 


1.98 ±0.65 


3.12 ± 1.10 


3.53 ±0.93 


Error in W + Acceptance (AA) 


Uncertainty 




Cut 1 


Cut 2 


Cut 3 


Missing NNLO 
Scale Dependence 




-0.24 ±0.90 
1.29 ±0.61 


-2.13 ± 1.18 
0.34 ± 0.71 


-2.53 ± 1.02 
0.20 ±0.73 


Total 




1.32 ±0.62 


2.15 ± 1.17 


2.54 ± 1.02 



Table 12: Summary of QCD uncertainties Act in the cross-sections and AA in the acceptances 
relative to the un-cut cross-section a tot . The three cuts are described in Table |l|. Missing NNLO is 
shown with a sign, because it has been calculated. 



Total Theoretical Uncertainty (%) 





w~ 


W+ 


Uncertainty 


Cross-Section Act 


Acceptance AA 


Cross-Section Act 


Acceptance AA 


Missing O(a) EWK 
Total QCD 
PDF 


3.91 ±0.57 
1.98 ±0.65 
3.31 


0.05 ± 0.54 
1.32 ±0.62 
2.22 


4.00 ±0.61 
1.66 ±0.69 
4.01 


0.09 ±0.59 
1.00 ±0.78 
2.28 


Total 


5.49 ± 0.47 


2.58 ±0.32 


5.90 ± 0.46 


2.49 ±0.31 



Table 13: Total theoretical uncertainty on the W production cross-section Act, and acceptances 
AA. 
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I W*->r v NLO Acceptance Cut 1 | 



| W*->r v Fractional NLO Acceptance Error Cut 1 | 




CTEQ6.5 MRST'06 (NNLO) CTEQ6.1 MBST'04 (NNLO)MRST'04 (NLO) MRST01 

PDFSet 



| W*->r v NLO Acceptance Cut 2 | 




CTEQ6.5 MRST'06 (NNLO) CTEQ6.1 MBST'04 (NNLO)MRST'04 (NLO) MBST'01 

PDFSet 



| W*->r v NLO Acceptance Cut 3 | 




CTEQ6.5 MRST'OS (NNLO) CTEQ6.1 MBST'04 (NNLO)MRST'04 (NLO) MBST'01 

PDFSet 




CTEQS.5 MRST'06 (NNLO) GTEQS.1 MRST'04 (NNLO )MBST'04 (NLO) MBST'01 

PDFSet 



(a) 



| W + ->r v Fractional NLO Acceptance Error Cut 2 | 




CTEQ6.5 MRST'06 (NNLO) GTEQ6.1 MRST'04 (NNLO)MBST'04 (NLO) MBST'01 

PDFSet 



(b) 



| W*->r v Fractional NLO Acceptance Error Cut 3 | 




CTEQ6.5 MRST'06 (NNLO) GTEQ6.1 MRST'04 (NNLO)MBST'04 (NLO) MBST'01 

PDFSet 



(c) 

Figure 16: Comparison ofW + — * l + v (£ = e or /i) acceptances A, with several recent PDF cal- 
culations for acceptance regions (a) Cut 1, (b) Cut 2, and (c) Cut 3, as defined in Table |l|. The 
left-hand plots show the total acceptance and the right hand plots show the fractional error on the 
acceptance. 
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I Leplon P T Cut2 | 




P T (GeV/c) 



| Lepton P T Cut3 | 




P T (QeV/c) 




(a) 



I Leplon P T Cut2 | 




(b) 



I Leplon P T Cut3 | 



(c) 



P T (GeV/c) 




P T (GeV/c) 



Figure 17: The W~ — > l~v cross-section a (£ — e or /u), as a function of the p T cut for acceptance 
regions (a) Cut 1, (b) Cut 2, and (c) Cut 3, as defined in Table [l[ For each acceptance region we 
fix the invariant mass and \r]\ cuts at their specified values, and vary only the p T cut. The figures 
on the right show the relative errors in the cross sections. 
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I Leplon P T Cut2 | 




| Lepton P T Cut3 | 



P T (QeV/c) 



35 40 45 




5 10 15 20 25 30 35 40 45 

P T (GeV/c) 



(a) 



| Leplon P T Cut2 | 




(b) 



I Leplon P T Cut3 | 




(c) 



P T (GeV/c) 



Figure 18: The W + — » l + v cross-section a (£ = e or fi), as a function of the p T cut for acceptance 
regions (a) Cut 1, (b) Cut 2, and (c) Cut 3, as defined in Table [l] For each acceptance region we 
fix the invariant mass and \r]\ cuts at their specified values, and vary only the p T cut. The figures 
on the right show the relative errors in the cross sections. 



|Lepton lr| I Cut2 | 




(c) 

Figure 19: The W~ — > Z 17 cross-section er (1= e or /i), as a function of the |ry| cut for acceptance 
regions (a) Cut 1, (b) Cut 2, and (c) Cut 3, as defined in Table [l[ For each acceptance region we 
fix the invariant mass and p T cuts at their specified values, and vary only the \rj\ cut. The figures 
on the right show the relative errors in the cross sections. 
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(c) 

Figure 20: The W + — > cross-section cr (£ — e or /i), as a function of the \r/\ cut for acceptance 
regions (a) Cut 1, (b) Cut 2, and (c) Cut 3, as defined in Table [l] For each acceptance region we 
fix the invariant mass and p T cuts at their specified values, and vary only the \rj\ cut. The figures 
on the right show the relative errors in the cross sections. 
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Figure 21: The W~ — > l~T7 acceptances (£ — e or /i) A, as a function of the p T cut for acceptance 
regions (a) Cut 1, (b) Cut 2, and (c) Cut 3, as defined in Table [l[ For each acceptance region we 
fix the missing energy and cuts at their specified values, and vary only the p T cut. The figures 
on the right show the relative errors in the acceptances. 
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| Leplon P T Cut2 | 
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Figure 22: The W + — » l + v acceptances (£ — e or /i) A, as a function of the p T cut for acceptance 
regions (a) Cut 1, (b) Cut 2, and (c) Cut 3, as defined in Table [l] For each acceptance region we 
fix the missing energy and \r]\ cuts at their specified values, and vary only the p T cut. The figures 
on the right show the relative errors in the acceptances. 
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(c) 

Figure 23: The W~ — > Z ~77 acceptances = e or /j,)A, as a function of the \r)\ cut for acceptance 
regions (a) Cut 1, (b) Cut 2, and (c) Cut 3, as defined in Table [l[ For each acceptance region we 
fix the missing energy and p T cuts at their specified values, and vary only the cut. The figures 
on the right show the relative errors in the acceptances. 
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(c) 

Figure 24: The W + — > l + v acceptances^ = e or fx) A, as a function of the |ry| cut for acceptance 
regions (a) Cut 1, (b) Cut 2, and (c) Cut 3, as defined in Table [l] For each acceptance region we 
fix the missing energy and p T cuts at their specified values, and vary only the \q\ cut. The figures 
on the right show the relative errors in the acceptances. 
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